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ABSTRACT 


The use of glass units as facade cladding in modern buildings has become a common practice. 
These units are generally formed as double or triple plates by dry air or another gas (argon, 
etc.) between the plates to provide heat insulation. This makes the engineering design under 
the wind and thermal loading difficult, since one must consider the force interactions between 
the plates due to the pressure change. In this study, the effect of the change in the geometrical 
properties of a doubly insulated rectangular glass unit on the interaction between the plates is 
investigated. Boyle-Mariotte law and geometric nonlinear plate finite element analysis are 
considered together. 


Keywords: Boyle-Mariotte's Law, Façade cladding, Finite Element Method, Geometrically 
nonlinear, Insulating glass units 


ÖZET 


Modern binalarda cephe giydirmesi olarak cam ünitelerinin kullanılması oldukça yaygın bir 
uygulama haline gelmiştir. Bu üniteler termal yalıtım sağlayabilmek üzere genellikle 
aralarında kuru hava veya bir başka gaz (argon vb.) bulunan çift veya üçlü plakalar şeklinde 
teşkil edilmektedirler. Bu durum, rüzgâr ve termal yükler altında yapılan plaka dizaynına 
ilişkin mühendislik hesaplamalarını zorlaştırmaktadır, zira plakalar arasındaki hacim ve 
basınç değişikliklerinden kaynaklanan kuvvet etkileşimlerinin deplasman ve gerilmeye olan 
etkilerinin dikkate alınması gerekmektedir. Bu çalışmada rüzgâr yükleri altında çift yalıtımlı 
dikdörtgen cam ünitelerine ait plakalar arasındaki etkileşimin ünitenin geometrik özellikleri 
ile değişimi parametrik olarak incelenmiştir. Hesaplamalarda Boyle-Mariotte yasası ve 
geometrik bakımdan doğrusal olmayan plak sonlu elemanlar analizleri birlikte ele alınmıştır. 


Anahtar kelimeler: Boyle-Mariotte Kanunu, Geometri bakımından doğrusal olmayan analiz, 
Giydirme cephe, Sonlu Elemanlar Yöntemi, Yalıtımlı cam ünitesi 


GİRİŞ 
Plaklar, bir boyutu diğer iki boyutu yanında oldukça küçük olan yapı elemanlarıdır ve cam 
ünitelerinin mekanik hesapları için kullanılması uygundur. Yalıtımlı cam (YC) üniteleri, 
birden fazla cam plağının termal yalıtım sağlayabilmek üzere aralarında kuru hava veya bir 


başka gaz olacak şekilde sızdırmaz teşkilinden oluşur. Yüksek binalar veya havaalanları gibi 
modem yapılarda genellikle yalıtımlı cam üniteleri tercih edilmektedir. 


YC ünitelerinin yapısal davranışını analiz etmek için literatürde birçok özgün çalışma 
mevcuttur. Chou ve ark. konu ile ilgili olarak gerçekleştirdikleri deneysel ve teorik 
çalışmalarında yalıtımlı cam ünitelerinin yapısal davranışları üzerine basitleştirici 
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formülasyonlar elde etmişlerdir (1 |. (2| de yalıtım cam ünitelerindeki yük paylaşımı, doğrusal 
ve doğrusal olmayan plak deplasmanları analizi gerçekleştirilerek ortaya konulmuştur. Kabuk 
eleman kullanılarak cam plakların doğrusal olmayan davranışını analiz etmek üzere sonlu 
eleman modelleri uygulanmıştır (31. Yapısal olmayan pencere cam plaklarının sismik 
davranışı incelenmiştir (4, 5J. Diğer araştırmacılar, ince cam plaklardaki yaklaşık yanal 
deplasmanları tahmin etmek amacıyla bir yöntem oluşturmuşturlar (6). Başka bir çalışmada 
ise birleştirilmiş değişken yükler altında yapısal camlar için bir doğrulanma formülü 
sunulmuştur (7). Çevresel yükler için yalıtım cam ünitelerinin yük paylaşımını belirlemek 
üzere karşılaştırmalı bir çalışma gerçekleştirilmiştir (8). 


Bu çalışmada, çift YC ünite plaklarının yayılı yük altındaki etkileşimlerini dikkate alan 
geometrik yaddoğrusal Sonlu Eleman analizleri gerçekleştirilmiş ve değişken yük ve yalıtım 
boşluğu için ünitenin mekanik davranışı parametrik olarak incelenmiştir. 


YALITIM ÜNİTELERİNİN SONLU ELEMAN FORMÜLASYONU 


Yükleme altında, YC üniteleri arasında hacim değişikliği nedeniyle kuvvet etkileşimi 
meydana gelmektedir. Bu etkileşim geometrik yaddoğrusal probleme ilaveten yük açısından 
da yaddoğrusal bir durumu getirmektedir. Problemin çözümü için Klasik Plak Teorisi (9| ve 
Boyle-Mariotte yasasının birlikte kullanılması gerekmektedir. Modellerde, elastik ve izotropik 
malzeme kabülü yapılmış olup, aynı zamanda yükleme esnasında YC ünitelerinin hem iç hem 
de dış kısımlarında sıcaklık değişikliği oluşmadığı ve ortam basıncının atmosferik basınca 
(Ppenv—-101.3 kPa) eşit olduğu kabülleri yapılmıştır. Şekil 1? de plak elemanına ait serbestlikler 
ve yükleme yönü gösterilmiştir. 
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Şekil 1. 4 nodlu kabuk elemanı 


Sırasıyla, yerel x , y ve z yönlerinde, u, v ve w yer değiştirme alanı olarak belirlenmiştir. 
Elemanın şekil değiştirme alanları ve bünye deklemleri Denklem (1) de verilmiştir; 
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burada 4, and P. sırasıyla uzama ve eğilme rijitlikleri olup, Denklem (2)'deki genel bünye 


denklemlerinden elde edilmişlerdir. 
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Sonlu eleman denklemleri, Virtuel İş Teoremi kullanarak Deklem (3) seklinde elde edilebilir; 
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Üst çizgili terimler yer değiştirme alanlarının varyasyonlarıdır. Şekil 2 “de, YC ünitesinin dış 
ve iç plakları ve seçilen yükleme yönü gösterilmiştir. 
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Şekil 2. İki tabakalı yalıtım cam ünitesi 


Dış ve iç plak denklemleri Denklem (3) deki tekil cam plaklar için kompakt formda aşağıdaki 
gibi yazılabilir, 


k,(W,)-g—d 
,(W,)g—dp (4) 
k,(W,) dp 
burada dp boşluk basıncındaki değişikliği (dp) temsilen oluşturulan sonlu eleman dış yük 


vektörüdür. İki plak arasındaki etkileşimi temsilen Boyle-Mariotte kanunu aşağıdaki şekilde 
elde edilebilir, 


k(dp) dp, KAV) p AV, <0 O) 


Burada V. ve p,,, sırasıyla, yükleme öncesi boşluk hacmi ve ortam basıncıdır. AK is boşluk 


hacimindeki değişikliği ifade etmektedir. Denklem (5)'e ilişkin Newton iterasyon şeması 
aşağıdaki şekilde yazılabilir. 


dk, 
ddp 


Analizler, Denklem (4) ve Denklem (6) birlikte kullanılarak iteratif olarak gerçekleştirilmiştir. 
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SAYISAL ÖRNEKLER 


Plak formulasyonun ve sınır koşullarının doğruluğunu test etmek üzere öncelikle boyutları 
2.0x2.0 (mxm) ve kalınlığı h-x6 mm olan tekil bir cam plak çözümü yapılarak sonuçlar Şekil 
3'te ANSYS ve pr-EN16612 ile karşılaştırmalı olarak verilmiştir. Sınır koşulu olarak: kenar 
nodlar plak düzlemine dik doğrultuda tutulu, orta nod ise plak sadece düzleminde tutulu 
şeklinde girilmiştir. Malzeme özellikleri için Hlastisite modülü E-70 GPa ve Poisson oranı 
v>0.22 olarak alınmıştır. 


(a) ANSYS (b) Bu çalışma (Denklem (3)) 


Sekil 3. 2x2 (mxm), h-6 (mm) tekil cam ünitesi Sonlu Eleman çözümü deplasman kontur 
grafiği (Eleman ağı 16x16). Yayılı yük g-1.0 kPa. Doğrusal analiz maksimum deplasman 
(W0): 48.78 mm. Geometrik yaddoğrusal analiz sonuçları (maksimum deplasman): 
ANSYS-19.13 mm, Buçalışma-18.43 mm, o pr-EN16612-19.45 mm. 


Parametrik incelemede hesap zamanını kısa tutmak üzere Şekil 3'te mevcut çözümde görülen 
“05 civarındaki hatanın kabul edilebilir sınırlar içinde olduğu değerlendirilmiştir. 


Tablo 1? de parametrik inceleme için seçilen YC üniteleri, yayılı yükler ve boşluk mesafeleri 
gösterilmektedir. 


Tablo 1. Parametrik inceleme için seçilen cam üniteleri (Sonlu Eleman ağı: 16x16) . 


Lx(m) Ly(m) hi/h2 (mm) g (kPa) c (cm) 
Ünite-1 1.5 1.5 6/6 0.6,1.0,1.4 o(0.5,1.0,1.5,2.0,3.0,4.0 
Ünite-? 3 1.5 12/6 0.6,1.0,1.4 o(0.5,1.0,1.5,2.0,3.0,4.0 


Analiz sonuçlar Şekil 4” deki grafiklerde verilmektedir. Bu grafiklerde “0” indisini barındıran 
deplasmanlar ilgili camın aynı yükleme altında tekil olarak çözümünden elde edilen 
maksimum deplasman değeridir. 
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Sekil 4. Tablo 1 de verilen YC ünitelerine ait analiz sonuçları (Wi: i numaralı camın 
maksimum deplasmanıdır). 


Grafiklerden de görüleceği üzere genel olarak c nin artmasıyla birlikte camlar arasındaki 
etkileşim azalma eğilimindedir. Şekil 4a da görülebileceği üzere eşit kalınlıklı plaklar 
dumunda yük paylasımı dış plakta 9650-60, iç plakta 9050-40 bandında gerçekleşmektedir. Bu 
durumda, mühendislik hesaplamaları açısında güvenli tarafta kalmak üzere eşit kalınlıklı 
plaklar için dış cama yayılı yükün *060' 1, iç cama “050” si gelecek şekilde tekil plak 
çözümleri ile hesaplama yapmanın uyun olacağı sonucuna varılabilir. Yine Şekil 4 (a) ve 4 (b) 
de görüldüğü üzere dış cam plak kalınlığı iki katına çıkarıldığında etkileşim büyük oranda 
azalmakta ve bu durumda oranlar dış plakta 9085-90, iç plakta 9012-22 bandına taşınmaktadır. 
Bu sonuç mühendislik hesaplamaları açısından 2/1 cam oranı için dış cama 9490, iç cama 422 
yük aktarımı ile tekil çözüm yapılabileceği şeklinde yorumlanabilir. Ara kalınlıklar için 
etkileşlimlerin h1/h2 oranına göre doğrusal değişeceği kabülüne göre hesap yapmak uygun 
olacaktır. 


DEĞERLENDİRME 


Bu çalışmada, YC ünitelerinin yayılı yükleme altında geometrik doğrusal olmayan analizleri, 
cam plaklar arasındaki etkileşim dikkate alınarak gerçekleştirilmiştir. Sonuçlar, plakların 
değişken kalınlıkları, ara boşluk mesafeleri ve basınç değeri için elde edilmiştir. Sayısal 
örneklerde yayılı yük ve ara boşluk mesafesinin ünite davranışına etkisi değerlendirilmiş olup 
tekil plak çözümlerinin ulaşılabilir olduğu durumlar için YC ünitelerinin mühendislik dizaynı 
açısından basitleştirici yaklaşık bir ön-dizayn kriteri önerilmiştir. Ele alınan problem, cam 
boyutları ve kalınlıklarının daha geniş data ile dahil edildiği bir incelemeye ihtiyaç 
duymaktadır. Bu inceleme sonucunda prEn16612 ye katkıda bulunacak formüllerin ve 
abakların geliştirlmesi bu çalışmanın devamını oluşturacaktır. 
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